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Abstract 

The i*2, Fq, R = Fl/Ft proton structure functions are derived in the QCD 
dipole picture. Assuming and renormalization-group factorization, we relate 
deep-inelastic proton scattering to deep-inelastic onium scattering. We get a three- 
parameter fit of the 1994 HI data in the low-x, moderate Q 2 range. The ratios 
Fg/F2 and R are predicted without further adjustment. The dipole picture of 
BFKL dynamics is shown to provide a relevant model for quantitatively describing 
the proton structure functions at HERA. The predictions for F<i and Fq are com- 
patible with the next-to-leading order DGLAP analysis, while R is expected to be 
significantly lower at very small x. 
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1. Introduction 



The purpose of the present paper is to show that the QCD dipole picture ^ |3], |j] 
provides a pertinent model for describing the proton structure functions at HERA in 
the low-x and moderate Q 2 range. Indeed the HERA experimental data on the proton 
structure function and particularly the most recent 1994 published data [||] give a good 
opportunity to check the validity of different QCD-inspired theoretical models. The QCD 
dipole model contains the physics of the BFKL Pomeron 0, i.e. the leading Inl/x 
resummation of the perturbative expansion of QCD. In order to apply this model to 
proton structure function phenomenology, we shall use the /c^-factorization properties |7] 
in the context of the dipole model. In this scheme we give predictions for the whole set of 
F 2 , Fq and R = F L /F T proton structure functions. F 2 is well reproduced with only three 
parameters and the ratios F G /F 2 and R are predicted without any further adjustment. 

To apply consistently the QCD dipole model to deep-inelastic scattering on a proton, 
we are led to introduce massive qq (onium) configurations inside the proton. On the one 
hand, we compute the coupling of gluons to the dipole states in an onium (see formulae 
([^[H])). On the other hand we exhibit factorization constraints on the non-perturbative 
probabilities of finding an onium into the proton (see formula (|19"D ). 

In section 2, we derive the expression for the onium structure function, including 
the determination of the gluon-dipole coupling using /^-factorization. in section 3, we 
use once more the fc T -factorization properties to get the predictions for the set of pro- 
ton structure functions Ft, Ft and Fq. They are obtained as inverse Mellin transforms 
which can be calculated by the steepest-descent method. In section 4, we adjust the 3 
free parameters of the resulting expressions with the 1994 data on F 2 , and give the pre- 
dictions for Fq and R without further adjustment. The summary of results is in section 5. 



2. Onium structure functions 



To get the structure function F 2 we rely on the /^-dependent factorization properties 
|7]] which are valid at high energy (small x). In a first step, we calculate the deep-inelastic 
cross-section a onium of a photon of virtuality Q 2 on an onium state. This onium state can 
be described by dipole configurations [fflj. The photon-onium cross-section reads 



g.onium 



J d 2 rdz<5> {0 \r,z)a(x,Q 2 ;r) (1) 



where $ < -°- ) (r, z) is the probability distribution of dipole configurations of transverse coor- 
dinate r (with modulus r) . 

In the /^-factorization scheme, one writes 

rQ' 2 _ r 1 dz 

Q 2 a(x, Q 2 - r) = d 2 k — &{x/z, k 2 /Q 2 ) F(z, (kr) 2 ) (2) 
J Jo z 

where a/Q 2 is the (7 g(k) — > q q) Born cross section for an off-shell gluon of transverse 
momentum k. F(z, (kr) 2 ) is the unintegrated gluon distribution of an onium state of size 
r and contains the physics of the BFKL pomeron. In the x-space, eq.(0) is a convolution 
. Using the inverse-Mellin transform 

1 rc+ioc 

f{x) = — duj x-»f(ou) (3) 

llTX Jc-ioo 



1 



2 



yields 

Q 2 a(x, Q 2 ; r) = ^- j 'duo J d 2 k &(u, ^)F(oo, (kr) 2 )x^ . (4) 
A second Mellin transform in the /c 2 -space, namely 

1 rc+ioo 

f( UJ ,(kr) 2 ) = — djf(u n )(kr)^ (5) 

ZlTT J c—ioo 

yields 

QV(x,Q 2 ;r) =J^J ^a(u,l-j)F(u n )(Qr) 2 ^. (6) 

Here we have used the Mellin transform of both a and F, and taken into account the 
relation 

J_^ ( ^-w +1 = w _ (1 _ 7)) . (7) 

The unintegrated gluon distribution for an onium of size r has been derived in the 
QCD dipole picture. One obtains |]J 



7T 7 uj - 2£*x( 7 ) 

where = 2*(1) - #(7) - *(1 - 7) with #(7) = and ^x{l) is the well-known 
BFKL || anomalous dimension, u (7) is, up to a factor, the Mellin transform of the 
gluon-dipole coupling, which expression in the QCD dipole framework we shall now give. 

Using once more the fc^-factorization in extracting a gluon of virtuality I from a dipole 
of transverse radius r, we find the following expression of v in the momentum space [fll 

v(rl) = — d9 (2 - e ilrcose - e"* cose ) = 1 - jMr) (9) 
Air Jo 

and correspondingly 

«(7) = / «(/r) (/O-^dC/r) = - r i + ^ - (10) 

A sketchy proof of the determination of v(lr) is obtained by using the eikonal coupling 
of the off-shell gluon to the dipole of transverse size r. We use the well-known light-cone 
coordinates as well as the light-cone gauge. When the emitted gluon is soft, the dominant 
component of its polarization is 

e> = f- (11) 

where e A is a two dimensional vector depending on the polarization index A. Given a quark 
(antiquark) of transverse coordinate x (x x ), the coupling has the form J^e MA ~ J+e-x 
with the corresponding eikonal current given by 

J+ = zgT a (e- d -^ - e-^° A, (12) 
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where a = g 2 /(47r). After averaging over color and polarization, this leads to the cross 
section integrand 

^-Tr\J + e x _\ 2 = ■^Ex=i 2 9 2 Tr(T a T a )(e- il --^ - e -fe)(e^ - e^)4(^) 2 

= g 2 C F (2 - e^o _ e-^io)^. (13) 

After integrating over the polar angle 6, and with the suitable normalization |7j, one gets 
formula @. Note that a similar expression is encountered in the elementary dipole-dipole 
cross section . 

Since F(u;, 7) (formula flD) exhibits a pole at u p = ^^xil)^ the lo integral yields 

Q 2 a(x,Q 2 ;r) = ^ [ ^( W = ^l- 7 )HM W )7 e ^xWMi). (14) 

Furthermore <r(u;, 1— 7) is weakly odependent [|7|], and thus 

a(u p , I-7) ~ <x(0, 1-7) = 4vr 2 a e . m ./i(7), (15) 

factorizing out the electromagnetic coupling. 

In order to average over the wave function of the onium state (see formula (|l|)), one 
defines 

< r 2 > 7 = I d 2 r (r 2 ydz $ (0) (r, z) = (M 2 ) -7 (16) 

where M 2 is a perturbative scale characterizing the average onium size. This means that 
M 2 is high enough to use the perturbative dipole-cascade mechanism [I]. The photon- 
onium structure function then reads 



7T 7 2Z7T \M 2 J V ' 7 

- /^fSV^,7), (17) 
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where F(x,j) is the inverse Mellin transform of the onium structure function. 



3. Proton structure functions 



Let us now deal with deep inelastic scattering on a proton target. Assuming the 
/cT-factorization properties to be valid for high-energy scattering off a proton target, we 
substitute in formula (|i7|) F(x, 7) by F(x, 7)^(7, M 2 ; Q 2 ) where to can be interpreted 
as the Mellin-transformed probability of finding an onium of transverse mass M 2 in the 
proton. Qq 2 is a typically non-perturbative proton scale. Thus, 

f— (*,Q 2 ;Ql) = ^/ ^M7)^ (^) 7 e^*"^» ro (7,M';^). (18) 

As previously mentionned, M 2 is an arbitrary perturbative scale. The overall result 
has to be independent of the factorization scale, provided it is in the perturbative region. 
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I IK' \ (UK 1. ! Ml ..1 / 

w has to match the (M 2 ) -7 dependence of the term (tL) 7 . One then writes 



Hence, assuming the renormalization group properties to be valid, the M 2 dependence of 

M 2 ■ 



fM 2 V 

w( 1 ,M 2 ] Q 2 )=w( 1 ) f-^J . (19) 

The argument leading to formula (|T9|) is similar to the one applied to the one- loop QCD 
evolution equation in ref. |T2[ . This yields the final result 

F^^Q 2 -Ql) = 2^/^(7)^(7) (|) 7 e^(7)in(i). (2Q) 

Note that F proton depends on the ratio Q 2 /Q 2 where Q 2 can be interpreted char- 
acteristic scale of the non perturbative hadron-onium structure function. 

Let us use the generic result (20) to derive the prediction for specific proton structure 
functions in the QCD dipole picture: 




( h 



n J 2m \Qc 



T 



v{l) 



h L -AIZ W ( 7 ) (21) 



V i 
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where Ft(l) is the structure function corresponding to transverse (longitudinal) photons 
and Fq the gluon structure function. The coefficient functions 



h T \ = a (r(i- 7 )r(i + 7)) 3 i / (i + 7 )(i - 1: 

h L J 3vr 7 r(2 - 2 7 )r(2 + 2 7 ) 1 - | 7 I 7(1 - 7) 



(22) 

are given in ref J7[. Note that the coupling constant a in formula ( ^2] ) is kept identical to 
the effective coupling of the BFKL kernel in formula (|8]). In the following we shall take 
a constant and not running, as in general considered in both BFKL and /^-factorization 
schemes. 

The integral in 7 is performed by the steepest descent method. Assuming that if (7) 
is smooth and regular near 7 = |, we can calculate the structure functions by expansion 
of the BFKL kernel near this value. The saddle point is at 

7 < = 2 (1 - <l " 1 % ) 



5JV P . t ,„ 4 , 1 



-1 



7C(3)ln- . (23) 



7T X, 

Formulae (^) show that the considered approximation is valid when 

a In A ~ hi -^-/ hi - < 1, 

that is the small x, moderate Q/Qq kinematical domain. 
We obtain 

F 2 = Ft + Fl = Ca i/2 e (ap-i)ini^ e -§m 2 % t 

Qo 



(24) 
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where 

ftp — 1 = (25) 

7T 

C, Qo and ftp are three unknown parameters which have to be determined by the fit. Once 
the F 2 fit performed, formula ( pO] ) allows one to predict Fq and R = Fi/Fp. Moreover 
R and Fq/F 2 are independent of the overall normalization C, i.e of the non-perturbative 
w(j c ). One gets 



F G 1 



F 2 /i T + h L 



= 37r 7c 1 - f 7 c r(2 - 2 7c )r(2 + 2 7c ) 

7=7c " a i + | 7c -|7c 2 (r(i-7c)r(i + 7 c)) 3 



(26) 



and 



p h L{ \ 7c(l - 7c) , 07 ^ 

= jfW = (i + 7o )(i-f)- (27) 

4. Comparison with experiment and predictions 

In order to test the accuracy of the F 2 parametrisation obtained in formula (|24]), a 
fit using the recently published data from the HI experiment |J has been performed. In 
order to remain in the domain of validity of the QCD dipole model, all the points with 
Q 2 < 150 GeV 2 were used in the fit (we checked that the result is not strongly dependent 
on the cut value). The parameters used for the fit are C, Qo, and ftp. The fit for F 2 is 
displayed in figure 1. The \ 2 is 101 for 130 points for parameter values: 

a P = 1.282 

Qo = 0.627 GeV. 

C = 0.077 (28) 
The value of the hard Pomeron intercept ap is in agreement with values found in other 



determinations using BFKL dynamics for describing F 2 at HERA fll3| . It corresponds 
to an effective coupling constant a ~ .11. As a matter of comparison, the input value 
corresponding to the HI QCD study is .12. As discussed in section 3 of this paper, Q 
and C are parameters of non-perturbative origin. The value of Qo, which corresponds to 
an average transverse size of .3 fermi, is in the range expected for deep inelastic scattering 
on a proton. C is connected to the non-perturbative input w(l/2). 

With only three parameters, the fit of the data is remarkably good in the low Q 2 range. 
When considering the whole set of HI data (see fig. 2), deviations from the fit are only 
observed in the high-Q 2 range (Q 2 > 200 GeV 2 ), for x > 10~ 2 . This is expected since the 
dipole model does not contain the valence contribution and does not describe the high Q 2 
QCD dynamics. By comparison with our previous results |10| based on the HI and Zeus 



1993 data ||, the quality of the present fit confirms the validity of our predictions, while 
taking into account the extended kinematical domain and the more precise data taken in 
1994. 

Using formula (|26"D , we get now a precise prediction of the gluon density. The value of 
a we used in formula (EBf) is the same as for F 2 , that is the one extracted from the fitted 



value of ftp by the help of eq.(^). The comparison with the published HI determination 
H is shown in figure 3. The grey band is the gluon value obtained by the HI collaboration 
from their DGLAP next-to-leading order QCD fit for two different values of Q 2 . We can 
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see that our prediction is quite satisfactory without any further adjustment for the HI 
values of Fq for x < 5.10 -2 (this parametrisation is not supposed to be valid for higher 
values of x). We notice that the fit is located in the range quoted by the HI collaboration, 
while all the gluon density values situated in the error band are equiprobable |5|, [|. In 
figure 3, we have also shown the result using only the first-order perturbative expansion 
of hr and hi in Eq. (p2|) . This corresponds to taking 



Fq 37T7 



F 2 0(1 + 7) 



(29) 



7=7c 



The comparison between Eqs. (|26| ) and ( p9|) exhibits the resummation effect of the In 1/x 
terms on the coefficient functions and h^. The result of formula (^) is shown as 
a dotted line in fig. 3. The HI determination clearly favours the resummed result. Our 
results show in particular that the gluon distribution does not allow to make a distinction 
between the DGLAP equations, used for the HI determination, and the BFKL dynamics 
expressed by the QCD dipole model. 

Let us comment the prediction of this model for the value of R = F^j Ft (formula |27|). 
We first notice that this prediction is independent of the non-perturbative normalisation 
C obtained in the fit of the structure function F 2 . The only non-perturbative parameter 
which enters the R prediction is Qq. The corresponding curve (full line) is displayed in 
figure 4. As before, it is useful to compare this prediction with the one-loop one (dotted 
line), namely R = 7 C . 

The measurement of R might be an opportunity to distinguish between the BFKL and 



DGLAP mechanisms. One expects the Altarelli-Martinelli [14| predictions for R to be 



larger [|TIJ than the limiting value 2/9 obtained in the dipole model. Moreover, the esti- 
mated values quoted in ref. ||, using the HI determination of the gluon structure function, 
are significantly higher than our predictions. It is known that this measurement is very 
difficult to achieve and the R values are very much correlated to the gluon densities when 
a QCD fit is performed. However, the fact that a sizeable difference for R is predicted in 
the same region where F 2 and Fq are well described is an incentive for further studies of R. 



5. Summary of results 



The main results of our analysis of deep-inelastic scattering can be summarized as 
follows 



i) Using (twice) the /^-factorization properties of high-energy hard interactions)]?! we 
were able to give a generic formula for the different structure functions of the proton, and 
by straightforward extension to any physical target, namely 

F T ^(X,Q 2 ;QD = & (^) 7 / i ( T )HM e ^x(7)ln(i) w Ta, ffet(7)) ( g 0) 

where h(j) is associated J7J to the (7 g — > q q) Born cross section (see (|ToT)), v (7) is the 
gluon-dipole coupling function (|i~0|) and w Target is the probability of finding an onium in 
the target. 
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ii) Using the steepest-descent method, we obtain a 3-parameter expression for F 2 which 
gives a fair description of all the 1994 HI data with Q 2 < 150 GeV 2 . 

hi) Without further adjustment, the model gives definite predictions for F G and R. 
One obtains: 



F G 1 



F 2 h T + h L 



7=7c "'T 



(31) 



7=7c 

where the expression of the saddle-point 7 C is given in formula (^3|). In particular, the 
resulting F G is fully compatible with the determination by Hl||. 

iv) The predictions of the QCD dipole model are similar to those obtained through 
next-to-leading order DGLAP evolution equations but for R, which differs from the pre- 
dictions at very small x using the Altarelli-Martinelli equation. 
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FIGURE CAPTIONS 



Figure 1 3-parameter fit of the HI proton structure function for Q 2 < 150 GeV 2 . The 
fit is constrained by the 1994 data from the HI experiment || (including the point at 2 
GeV 2 not shown on the figure for convenience). 

Figure 2 Comparison of the fit with all 1994 HI data. The validity of the fit extends 
beyond the domain included in the chi-squared (the highest Q 2 point at 5000 GeV 2 is not 
displayed for convenience). We note a discrepancy at high x, and high Q 2 , as expected 
outside the validity of the model. 

Figure 3 Prediction for the gluon density Fq(x, Q 2 ). We compare our prediction (full 
line) (cf. text) with the gluon density obtained by the HI collaboration |J (grey band) 
from their DGLAP next-to-leading order QCD fit. In order to exhibit the resummation 
effect, we also display the one loop approximation (dotted line), cf. formula (|29|). 

Figure 4 Prediction for R. The full line describes our prediction, and is compared to 
the one-loop approximation (dotted line). 
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